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§ ; Abstract 

The method has been developed to calculate effects of polarization not only 
for a atomic core in a field of valent electron, but also polarization of atom as a 
whole in the electron-hole formalism. A secondary quantized density matrix for 
many-electron system was used to find the Green function of a quasiparticle and 
its effective mass due to many-particle effects. 
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1 Introduction 

m 

As is known p], for odd number of electrons the Hartree - Fock method of self-consistent 
field gives the equation, calculation with which help is connected to a problem of diago- 
nalizing a matrix Lagrange multiplier. It results to possible noncommutativity of a Fock 
operator with the operator density matrix (a projector on subspace of orbitals for electron 
in a atomic core) and, hence, to absence of an interpretation of the Fock operator as a 
Hamiltonian for single-particle state [2], [3], [I]. The method of pseudo-potential allows to 
qh' interpret diagonal elements of the Lagrange multiplier in a Hartree - Fock equation as the 
energy of single-particle state. However, a pseudo-valent orbital constructed by the method 
appears so compressed one, that results to very underestimated lengths of chemical bonds in 
comparison with experiments [5], [B]. It, seemingly, is stipulated by the incorrect description 
of polarization effects in an atom as in the method of pseudo-potential one uses approxima- 
tion of spherically symmetric unexcited atomic core ("frozen" core). Therefore, at the last 
time in the "ab initio " calculations fulfilled the pseudo-potential is utilized only within the 
limits of core, supplementing with the phenomenological potential describing polarization of 
the core in the classical way [7j, [8]. Although, for example, for halogen dimers a relativis- 
tic calculation taking into account the dipole polarizability of atomic cores gives lengths of 
chemical bounds matching satisfactory with experiment data, the found binding energy is 
underestimated with respect to its real value [9], [TU] . [TT] . Obviously, it is due to the fact 
that the correct quantum description of polarization of an atom is possible, if one assumes 
the existence of quasiparticle exitations for the spherically symmetric core. 

The goal of this paper is to develop a self-consistent method allowing to calculate effects 
of core polarization in a field of valent electron within the framework of secondary quantized 
electron-" hole" formalizm. 



2 Wave function of many-electron system 

Let an arbitrary function ip to depend on coordinates r*j, i,...,n of n electrons: tp = 
^iXit • • • j ^n)- To be of a wave function for a system from n electrons the function 
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i^i^an^azi ■ ■ ■ !*«„) has to satisfy the Pauli principle, i. e. to be antisymmetric over space 
coordinates. Here the index oij, i — 1, . . . , n runs {1,2,..., n} so that on ^ oij for i ^ j. It 
can be achieved if one represents this function as 

*P(r ai ,r a2 ,...,r a J = e(P aia2 ___ a Jip(n,f 2 , . . . ,f n ), (1) 

where P ai a 2 ...ot n is a permutation: 

PaiQ2 - a " = ( ai a 2 . . . a n ) (2) 

which converts 1 into «i, 2 into a 2 and so on. A symbol e(P) denotes a number equal to +1 
if the permutation is even, and a number equal to —1 if the permutation is odd. 
It is easy to prove that the functions p]) satisfy the following equality: 

<*!,?<*„..., raj = 0. (3) 

To do it we divide up the left side of the expression into functions so that one describes 
the spin configuration as {ft ••• TT I III ••• !}> an d another one describes the spin con- 
figurations: {TT • • • Tl I Til • • • I}, {TT • ■ ■ Tl I ITI ■ ■ ■ I}, {TT • • • Tl I IIT • • • I}, 
{TT • • • Tl I III • • • IT}) with the sign minus "— ". Graphically this partitioning is shown in 
fig. 1 and is written mathematically as 

ip(ri, . . . ,fk-i,rk\rk+i,r k +2 ■ ■ ■ ,rj = ip(ri, • • . ,r k _i,r k+1 \rkr k+2 . . . ,f n ) + . . . 

+i>(ri, • • • , r fc _i, r k+ i\f k+1 , r k+l _ u f k , r k+l+1 , . . . , f n ) + . . . 

+^(ri, . . . ,r fc _i,r n |f fc+ i, . . . ,f fc _ n ,f fc ). (4) 



{ tttt i mil} = { mi I mil } + { mi I ij^} + {mi 1 1^} + --- + {mil mm) 

k n — k k n — k k n — k k n — k k n — k 

Fig. 1. Graphics represents the property of cyclic symmetry for a wave function of electron. 

But the wave function situated in the right side of symbolical expression in fig. 1, describes 
a configuration obtained by a cyclic permutation of electrons from the configuration at the 
left in fig. 1. Thus we have proved, that the electron function is symmetric with respect to 
the cyclic permutation and a mathematical notation of this property of cyclic symmetry is 
the expression (jlj). 

The set of functions (JTJ) is a basic set to construct the wave function of a many-electron 
system. Slater determinants possess the properties such as functions from the introduced 
basic set. 

Further we use a method of secondary quantization to establish a Hartree - Fock equation 
describing a single-electron state in the self-consistent field, within the framework of electron- 
hole formalism. 



3 Electron-hole formalism 

Hamiltonian H for the many-electron system reads: 



E 

i=i 



2m 



V 2 (fi) + V(fj) 



E 



TV — r~\ 

i>j=l 11 J' 



i>j=i 



i=l 
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where V 2 is the Laplacian, V(f£) is the potential energy of i-th electron in a field of atomic 
nucleus or nuclei, V(fi — r}) is the Coulomb interaction potential. 

Let us go into representation of secondary quantization where a single-particle state is 
given by creation operators ^(xi) and annihilation operators ip(xi) of i-th Fermi particle 
with generalized coordinates X{ = {r*j,tj,<7j} being its radius-vector r i; time £, and spin <7j. 
These operators satisfy the commutation relations |12j : 

^(x')^j\x) + ifl (x)ip(x') = 5(x — x') (6) 
$(a/)$(x) + ${x)$(x') = 0. (7) 

Now one can introduce an operator of "hole" creation in the following way. Since the wave 
function of the system can be written in the form 

(*»> ■ ■ ■ > f k+2, r k+1 )^(f k , r fc _ 1; . . . , ri)|0) 
= [ini{r n ) + i>(n+i)i(r n +i)$l n - k+1)l (r n +i, • • • , r k+2 , r fc+1 )^J T (f fc , r fc _ 1; . . . , fi)|0) 

= fc(r n ) +'0( n+ i) T (f n+ i)]|'01,...,'0n+l), (8) 

the operator [^ n j(r n ) + ^( n+ i)|(r„ + i)] is a creation operator of a "hole" by n-th electron. 

Here ^ n Ar n . . . , f k+2} r k+ i)ijjl^(r k} r k -i, . . . , rx) is a secondary quantized wave function of the 
system describing the configuration from k electrons with spin "up" and n — k electrons with 
spin "down" , n = 2k + 1, i.e. there is one unpaired electron; |0) is a vacuum state. The 
wave function (jSJ) describes systems with outer unpaired electron or with a "hole", having 
spin "up" | by outer electron. 

Let us assume that core polarization happens at the moment t. Then in this moment 
t the secondary quantized wave function of system can be obtained as a result of cyclic 
permutation p( cycl >(t): 

P {CyCl \m nl (f n ) + f (n+1)T (f n+ i)]|^l, • • . , Vn+l). (9) 

Since the operators p( c ^ c/ )(t) and [ip n i(r n ) + ^(n+i)t(^n+i)] commute with each other then 
according to the definition of permutation operator which is shown in fig. 1 the expression 
(HOP can be rewritten as 

P {cyd \t)$ nl (r n ) + ^ (n+1)T (f n+ x)]|Vi, • • • , Vn+l) 
= hMr n ) + ^ n+1)] {r n+l )}P^ d \t)\^ Vw> 



Cnm(t)^ m l(r m , t) + ^ C„ m (t)Vv + i )T (f m+1 ,t) 



,m=l m=k+l 



C nm {t)^{x m )P^ d Ml, ^n+l), (10) 



m=l 



where m-th term of the sum describes a "hole" which is created at m-th core electron after 
the cyclic permutation Pm yd ^ in the sense of the expression (j3J), a matrix |c nm (t)| converts 
the creation operator of the "hole" by n-th electron into the creation operator of the "hole" 
by m-th electron. 

According to the property of cyclic symmetry (j3J) one can formally write the expression 

n 

fc(r„) + ^ n+1 ) T (f„ +1 )] = c nm (t)ij(x m )P£ ycl) . (11) 

m=l 
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It follows from here that the sum J2m=i c nm(t)4>(%m) has to satisfy the same quantum equa- 
tions of motion as for [^VlX^) + i)(n+i)t( r *n+i)]- Heisenberg equation of motion for the "hole" 
creation operator [VVlX^n) + ' l P(n+i)i(r n+ x)] reads: 

(12) 
(13) 

(14) 
(15) 



7^*$»»i( f n) + ?(n+l)T(*n+l)] = [fe(*n) + f(n+l)T (*n+l)] > ^ 

K=[it l f H(xi)dfi+ it, f f Vixux^dndfA 
\i=i J i>j=i J J ) 



1 



V{xi,Xj) = ^(x j )i/; i (x i )V(f i - FjWixiMxj) 

where [A, B] is a commutator of operators A and B, is the Hamiltonian operator in the 
secondary quantization formalism which is obtained by a secondary quantization procedure 
from the operator (jSj). 

Substituting the expression (ITT1) into the equation of motion ffl2l one gets 



hd_ 



n n ( „ 

J2 C nm(t)Hx m )P£ Vd) = Yl C nm(t)\ J df. 



m=l 
n 



m,i=l 



i=i 



il)\x i )H{f i )i{x i ) + ^ / -^■V t (2 ; i)V' t (^)V(fi - fj^x^fixj) 



i=i 



ffijx^H^fjxi) + ^ / -^^ t (a; i )^ t (^)F(f i -f i )'0(s i )-0(a; i ) 



^(x m )|p^ c/ ) (16) 
for j < z 

for the "hole" creation operator Ylm=i c nm(t)4>(xm) which depends on configuration. Using 
the commutation rules (j()J [7J for quantized fermion fields one transforms eq. fll6p to the form: 

h d 



xi/)(x m ) 



It lb I n 

m=l m,i=l L 



fe ra )Jp^ (17) 
for j < i 



where — x m ) is the Dirac 5-function. Differentiating over time at the left side and 
integrating over the Dirac 5-function at the right side of eq. (JIT]) one gets finally 

hd In Cn m (t) ~\ 



dt 



el) $(x m )P<w*> = \ H{r m )i{x m ) + J d ? 



H{r m )${x m )-J2 J d?3\xi)V(n-? m )$(x m )$(x^ 



(18) 
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where I is the unity operator, (Tis the operator of "hole" energy as 

$(x m ) =i}(f m ,a m )exp(-ielt/H); 

and the right side of eq. (fl8l) is rewritten taking into account matrix multiplication rules. 

Now we can find the equations, describing single-particle state neglecting correlations in 
movement of electrons relative to each other. To examine the configuration shown in fig. 1, 
we shall assume that electrons with spins "up" move independent on electrons with spins 
"down". In other words their movement are not correlated. Therefore the wave function of 
such configuration is factorized in the following way: 

${ n - k)i {r n , . . . jffc+^rfc+i^^ffcjr*-!, . . . ,n)\0) = |-0i, Vn) 

= 9[ n -k)i(fn, ■ ■ • ,^ +2 ,f fc+1 )|0 i)$J T (ffc,f fc _i,...,fi)|0 t). (19) 

From here it follows an expansion for the vacuum state |0) 

|0) = |0|)|0T) = \0,a m )\0,-<Ti), (20) 

which means that |0) consists of non-occupation states with spins "down" |0 J.) and with 
spins "up" |0 |). 

Hermitian conjugation of eq. fllSH has the form: 



= P^ ^H{r m )^{x m ) d^(x m )V(fi - f^ix^ixi) j . (21) 

Acting with the hermitian conjugate equation (j2ip on the found non-occupied states \0a m ) 
in the vacuum state (1201 . one gets the following equation 

= Pj*& (H(f m )^Jf m ) -f^J dr&jr m )ft ai (n)V(? t - ?m)$-*M j loOlo, 

(22) 

According to the expression (j3J) the operation of permutation Pm yc ^ entered into eq. (1221) is 
written in the explicit form as 

P^fojPm)ft* t W$^(?i) = ^(^)^ CTl (r i )^ m (rl) - e ra(fm )^t_ CTi(f . ) ^ (f . ) . (23) 
Substitution of the explicit expression for Pi cyc ° (TJHD into eq. (j22J gives the following equa- 



tion: 



ih 



dl nc nm (t) 
dt 



x (® i (»Vn)V , (fi - TVn)^I ffi (n)& m (n) - fijO^ ^ *m)^l ffi (fi)^ . i (ri)) 

x|0(7 m )|0, -<7j). 

Multiplying on the left side of eq. (I24|) by the vector (0, — (7j| one gets 
ift ^U^t ro ( fm )| 0(Tm ) _ (0,-^let 110,-^)^(0100 = H(f m )^ am (f m )\Oa m ) 

n „ 

-^2 ^v5i(Oi°o^(n - o(0) -^1^-^(^)^(^)10,-0-4) 



(24) 



(25) 



as (0, — cTj|0, — <7j) = 1. If one introduces the following designation: tpl (fk)\0<Jj) = ipj(xk) 
and represents the unity operator in the explicit form: 

n n 

1 = E^(^)i ^)( ^fe(^) = 

then eq. fl25l) can be rewritten as 

ih 9 ^ nC ™ t -0m (O - E e + Pj^ m (x m ) = H(r m )ip m (x m ) 

n „ 

+ / rifj (V>m(0^(^ - f m )ip*(xi)ipi(xi) - ipi(x m )V(fi - r m )tfj* (xi)ip m (xi)) . (26) 

i=l ^ 

Eq. ([2"B"1) taking at initial time i = is the equation describing single-particle state ip m (x m ): 



H(f m ) + V sc (x m )-t x (x r 



3=1 



(27) 



where the differentiation over time t taking at initial time t = is designated as e m (0): 

_dlnc nm (t) 



e m (0) = - %h- 



dt 



i=0 



V sc and S x are the Coulomb and exchange interactions, respectively: 

n „ 

V sc (xi)ifj n {xi) = ^2 ^m{ x i) v {\ri ~ rj\)^ m {xj) dr^ n {xi) : 

m=l 
n „ 

t, X (Xj)lj) n (Xi) = ^ / ^(x-,>(|fi - fjD^n^j) dfj1p m (Xi). 



m=l 



(28) 
(29) 
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Physical meaning of the operators (|28| 129]) becomes evident if one rewrites them in terms of 
spinless electronic density p(r, r 1 ) and supposes that the interaction v is the Coulomb one: 

= EfcJ^tCfllMO, u = e 2 /|r-r-|. (30) 

From here it follows that the operator V sc represents electrostatic interaction of one electron 
with the electron density produced by remaining n — 1 electrons and electrostatic self-action 
(s.a.): 

(n-l)/2 . 

V sc {xi)ifj n (xi) = ^ / ^m{rj)v{\ri - fj\)if> m (fj) dr^ n {fi) + s.a. 

rff /ifr'^V n^i) + s.a. (31) 



m=l 



l r j | 

Analogously one obtains that the operator E x gives quantum exchange with exchange self- 
action (s.a.): 

n—l „ 

t X (Xi)lf} n {Xi) = ^ dr j' l l J *m(rj^ (T j) v (\^i -^j\)^n{rj,(Tj)^ m {ri,(Ti) + s.a. 

m=l ^ 

n-1 

""(O-i-CTi] 



m=l 

+V4(^i, -Vj)i>m{rj, (7j)8{aj - ai))v{\fi - r}|)^ n (r}, aj) + s.a. 

d f e P( r i>ri) ^, ) + B a (32) 
|r< - r,-| 

Since the operators V rsc and into the expression (1271) are subtracted from each other then 
the self-acting terms (s.a.) entering into it are mutually cancelled. 

Let us suppose that there exists a representation in which the energy operator of a 
"hole" ?t entering into eq. (|27|) is diagonalized = e(ki)I, where I is an unitary matrix. 
The replacing f m — > r*j and the taking into account of the diagonalization condition into 
eq. (12?]) allow to describe the core polarization as a quasiparticle exitation with the energy 
e(ki) which stationary state is determined as 

H(fi) + V sc (k iXi ) - t x {k lX ^ = U m (0) + fy p i j MkiXi), (33) 

e+ = e(ki)I. (34) 

Approximation of a "frozen" atom 

Let us call by approximation of a "frozen" atom the calculations with assumption that 
almost all time a "hole" is near m-th electron which is located in a point with radius- vector 
r*j. Since the energy operator eo/"hole" by m-th electron in the point r*j can be of a work 
only which has to be perfomed to shift an electron occupying the "hole" near m-th electron 
to j-th non-occupied orbital of atom "frozen" in the moment, then e is equal to: 

^ = -^-4 m) )- (35) 
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Substituting the expression (1351) into eq. (1331) one gets 

~H{n) + V sc (x t ) - Z*{xi)] ^ m {xi) = ^e m (0) - f^e^ - e^P^j ^ m (x,). (36) 

The approximation of valent electron tp v or the approximation of a "frozen" atomic core 
appears if there exists an outer valent v-th electron located in a point with a radius - vector 
f v such that it always is outside limits of atomic core: r v > R c , where \R C \ is a radius of the 
"frozen" core. Existence of outer electron means that one of electrons can not be placed 
inside of the core. It is possible if the core has no "holes". Since in this approximation the 
core has no "holes" then eq. ( l36l) can be rewritten as 



HtrA+VotxA-TFtxi) 



Ipm(Xi) = ( e m (0) - (1 - Sic) _ 6 j) P 3 J ^rn(Xi), (37) 



the valent orbital is the orbital determined by the expression 

1p v = (1 - ^ Pc)i>m(Xi), i ^ C. (3^ 



Here index "c" labels core electrons. It is easy see that eq. fl37|) is the equation for Phillips - 
Kleinman pseudo-potential V PK [2], [3], [1] in the case of wave function depending on both 
coordinates and spin variables of electron. Therefore eq. (I3"3"|) and eq. (I3"B"1) are generalization 
of eq. (I3TI) and with their help one can describe many-electron systems beyond the framework 
where the method of pseudo-potential is valid. 

4 Secondary quantized density matrix 

We rewrite (I3"3"|) in the representation of the Dirac ket (bra) - vectors: 

N 



h{k)\n-k)+ J S{k-k') dk'{\n;k){m;k'\v{kk')\m;k') 

+(\n; k')5 nm )(m; k'\v(kk')(\m; k)S mn )) 

N 

-22 \m; k)(m; k'\v{kk')\n; k')S{k - k') dk' = |n; fc)(e n (0) + e n {k)) 

m=l J 



(39) 



where ki = {ki,<Ji}, h{k) is a momentum representation of the non-perturbed hamiltonian, 
v(kk') = J dfdf'\r)(k ■ r\v(\f— r*|)|fc' • f f )(f f \ is a momentum representation of the Coulomb 
interaction operator, S(k — k') is the Dirac 5-function manifesting the presence of the law of 
conservation of momentum. 

Let us introduce projective operators 



p^=\m-k'){n-k\- (40) 

and express eq. (I3"9"|) through these operators p™£. For this purpose, Eq. (I3"9~|) is multiplied 
from the right by bra- vector (n;k\. Then, additional summating over n and integrating over 
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dk lead to the equation: 

N „ N 



[ dk h{k)\n; k)(n; k\ + J^ [ / S ( k ~ k ') dk dk ' x 

n=l n=l ^ J 

k)^2(m;k'\v(k,k')\m;k')(n;k\+ j j 5{k - k') dk dk' x 

m=l ^ ^ 



JV \ / N 




^ |n; k'}5 nm (m; k'\ ) I ^ |m; k)5 mn (n; k\ 

m=l / \n=l 

N 

K k )( m ; k'\v{kk') \ n ; k ')( n ; k H k ' ~ k)d{-k)d{-k') 

N 

k\n; k)(e n (0) + e n (k)). 



m=l n=l 

N 



n=l 



(41) 



We see that the first and second terms in the left side of the equation (141]) are traces of a 
matrix representation of operators ph and p*pv , and the third and fourth terms in the left 
side of equation (14"I!) are mutually canceled. Hence, it means that using normalization of the 
function |n; k): j dk{n; k\n; k) = 1, we obtain the following equation: 

r N 

Spph + Spp*pv = e n (0)N + dk^2(n;k\n;k)e n (k). (42) 

J n=l 



Using properties of the projective operators p™$\ (p™$)* = p™$ and (p^) 2 = p™$, one can 



transform eq. ff42|) to the form: 

N 



Spp(h + v)= e n (0)N + / dk^2(n; k\n; k)e n {k) = e n (0)N + e. 

J n=l 



(43) 



Let us elucidate physical sense of introduced projective operators p™$. The energy of a 
quasiparticle e is in the right side of Eq. (H3|) up to the constant e n (0)iV. From here it follows 
that the operator p™Q allows to calculate the energy e of quasiparticle excitations. It means 
that the expression (j43p is nothing else but a procedure of average over density matrix. 
Since the averaging with the help of the operator p^f yields the energy e of quasiparticle, 
this operator is a secondary quantized density matrix. 

Thus, it is proved that eq. fl33l) can be considered as the equation describing the state 
of quasiparticle and determining its energy with accuracy up to the constant e n (0)iV. Since 
for description of single-electron state it is necessary to take into account the presence of 
"hole", then the quasiparticle state describes the electron- hole pair. 



5 Green function for single-particle state 

It follows also from Eq. fT45|) that the quantity e n (ki) can be interpreted as an eigenvalue 
of the Hamiltonian for a quasiparticle excitation without taking into account interaction 
of quasiparticles. Therefore, the equation (I43p . written in the formalism of the density 
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matrix Pni'-kk' = Pw can ^ e rewritten in the formalism of wave functions in coordinate 
representation and in the limit of large N, N — > oo in the following way: 

- (h + E* + V SC )) J^Pnl^ = ^J^n(0))N5 rr , , (44) 
' n 

where property J dk(n;kr\n;kr') = 5 rr > has been used; 5 rr > is a delta symbol. Since the 
energy e n (0) of the bound one-electron state is negative: e n (0) < 0, the right side of Eq. (jHJ) 
represents itself the Dirac (^-function 8(r — r'). It allows to write Eq. (lUi) as: 

4-h HF Wt%ry = 6(r-r'), (45) 

' n 

where h HF = (h + if - + V sc ), if - = — S x . Eq. (j4"5l) is the equation for the Green function. It 
means that in the secondary quantized representation the operator 

Gf\n';r,r') = J2p ( Xr' (46) 

n 

possesses properties of non-perturbed Green function. 

So, the quasiparticle excitation determined by the Hamiltonian h HF can be considered 
as a free particle whose equation of motion is the equation (|45p . 

In many-body problem, in particular, in calculations of energy-band of crystal structure 
a contribution given by the interaction of an electromagnetic field with matter plays an 
essential role. To take into account many-particle effects due to correlated motion of an 
electron we should describe the system by self-consistent solutions of the non- stationary 
equation 

%9 ^r = (47) 

where H is the Schrodinger hamiltonian in a non-relativistic case or a Dirac hamiltonian in 
relativistic one. 

We have proved that for the secondary quantized representation the operator p looks 
like p = and possesses the properties of the Green function G\. Therefore the sum 

G\ over n from elements of the matrix p^, for the secondary quantized density matrix p 
describing an interacting particle satisfies the Dyson equation in nonrelativistic case or to 
the Schwinger - Dyson equation in relativistic case: 

Gi(l;2) = GS 0) (1;2) + J d3 dA Gf\l; 3)11(3, 4)Gi(4; 2) (48) 

where G^\l; 2) is the free Green function, E(3, 4) is the self-energy operator: E = if + S c , 
S c is correlation interactions, representing itself a part of the self-energy which describes the 
many-particle effects. Here numerical labels for the arguments are used: {ri,t\} — X\ = 1, 
etc. Acting on Eq. f|48|) by the operator — h HF and using the equation of motion for the 
free particle (j4"5l) we get the equation for the perturbed Green function as 

'Ft ~ kHF{ri) X ' 2) ~ / d3tC ( n '> X ' 3 ) Gl ^' ; 3 ' 2 ) = (" e n(0))iV5 rir2 . (49) 
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Rewriting Eq. (j4"51) in the formalism of wave functions one gets 

- h HF (n) Mhn) - I df 2 ± c (n; 1, 2)V>„(fcir 2 ) = (-e n (0))^ n (hn). 



(50) 



Since the expression: i^f- = e n (ki) takes place, then Eq. f[5"0~j) yields the Hartree - Fock 



equation taking into account of interacting quasiparticles 

h HF (ri)Mhri) + [ df 2 t c (n;l, 2)^ n {k X T 2 ) = (e„(0) + e n (fc 1 ))^„(A;ir 1 ). 



(51) 



Let us define a mass operator AM as: 



dr 2 t\n- l,2)^„(fcir 2 ). 



(52) 



Within the framework of the concept of quasiparticle excitations it is possible to represent 
the operator AM in the diagonal form: 



•(AM n (0) + AM„(A; i ))V„(fc i r i 



(53) 



the eigenvalue of mass operator possesses the property: 

AM n (ki) = AM n (—ki). Here AM„(0) is an eigenvalue of mass operator AM in the limit 

k -> 0. 

From here it follows physical sense of AM. It determines an effective mass of the quasi- 
particle and an efficient bottom of the energy band: 



h M * (n^hn) = (e n (0) + i n (h))Mhri) = 
[(e n (0) + AM n (0)) + {e n {kx) + AM n (fci))] ^ n {k x n). 



(54) 



Extremum of energy band Extr E n (ki) in the presence of interaction between quasiparticles 
is determined by the following expression: 



ExtTEM = (AM„(0) + e n (0))iV. 



(55) 



Let us consider the Green function normalized per unit volume V — 1 so that average 
energy in V is equal to the energy of single-particle state and N = 1. If — e n (0) — > oo then 
Eq. (jl9]) describes propagation of single particle and should be rewritten as 



Gi(n';l,2) 



rf3S c (n';l,3)G 1 (n , ;3,2) = 
= (-e n (0))5 ri7 



(56) 



Then according to the definition of the Green functions we have the following expression for 
the energy e n (0): 



e n (0) = C-a r 



C — > oo; 



(57) 



where a n is a finite quantity. Hence, since the energy is counted off from an arbitrary value, 
Eq. fl55|) yields the following expression for the reference points e(0)^ of quasiparticle energy 
and antiquasiparticle energy 



e(0)± = ±a n = (ExtiEih) t AM n (0)J /2. 



(58) 
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Here one took into account that N = 1; an extremum of zone is redefined as ExtrE(ki) = 
ExtvE n (kx) — C, the sign {±} in the left-hand side denotes cases of quasiparticles and 
antiquasiparticles, respectively; and the energy of particles in the pair is counted off from 
zero level. From the expression ( 1581) one gets that a n is the energy which is required to create 
a pair from quasiparticle and antiquasiparticle when k% = because 

a n = (e(0)+ - c(0)~)/2. (59) 

Because of an additional term e n (0) in the right-hand side of Eq. fl54l) we, generally 
speaking, cannot examine the left-hand side as a Hamiltonian operator of quasiparticles 
system acting on a corresponding wave function and as a consequence, can not construct a 
basis set of single-particle states of the problem. However, further we show, that h HF is a 
Hamiltonian of an electron - hole pair. 

Non-relativistic case 

In non-relativistic limit one can examine quantum systems which are characterized by a 
small value of AM n (0): 

AM n (0) -> 0. (60) 

It means that weak many-particle effects occur and, accordingly, we can speak about a 
"light" electron. The equality fl58|) occurs under condition of f l60l only in the case if a n — 0. 
From here it follows, that the energy a n of a pair is equal to zero. In other words, the energy 
is not expended to create an electron - hole pair. 

Substituting Eqs. flBTl) . (I6"0"j) into flB^|) and taking into account the condition a n = one 
gets the Schrodinger equation in the form 

h HF ( ri )ij n ( ri ) = inMn), (61) 

which describes the quasiparticle - antiquasiparticle pair (a non-relativistic electron - hole 
pair). Here e n = e n — C . 

Since the energy a n , expended on creation of a pair, equals to zero we have proved that 
the variable e n can be understood as the energy of an electron - hole pair. Therefore, Eq. (IBTj) 
has a group of dynamic symmetry, which algebra is so(3) x so(3) ~so(4) if to neglect exchange 
interaction. As is known, a nonrelativistic hydrogen-like atom possesses such a symmetry. 
Hence, we have proved that to calculate quasiparticle states in non-relativistic case it is 
possible to use a basis set of states of a nonrelativistic hydrogen-like atom. 

However, for a heavy electron AM n (0) > 1 according to the formula (loT?]) we always have 

a n = -AM n (0)/2 (62) 

and, hence, there is no an equation such as the Schrodinger one for its description. From 
here we conclude that the heavy electron can not be examined in nonrelativistic limit. 

Relativistic case 

Let us generalize the proposed approach to relativistic case. To do it we substitute 
Eqs. (15Tj) and (16"2"j) into (151| and let n tends to n — ► oo: 

h HF ( ri )ij n (k iri ) = ( Aii y 0) + i n {h)j Mhn), n -> oo. (63) 

Then, one can assume that the operator -§[ — h HF in Eq. ( 1631) is a quasirelativistic hamiltonian 
written in the implicit form in the Hartee - Fock approximation. 

From consideration carried out above it follows that the desired relativistic equation 
of motion should describe a charged composite system from a pair of particles and have 
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the dynamic symmetry S0(4). A spin of a given quantum system should be equal 1 as 
motion of a hole is a motion of an electron in many-particle positively charged matrix. 
In [H] the equation of motion of a relativistic charged vector boson has been found and 
it was shown, that it describes a relativistic hydrogen-like atom. The relativistic charged 
vector-boson appears as a composite system with a corresponding spectrum of masses and 
in quasirelativistic limit n — > oo its energy E\ is determined by the expression: 

m rwy 2 rrry 4 / 4 3\ rwy 6 / 3 8 4 \ g 

1 ~ 2" ~ 2^ ~ 8^VM ~ ^\y~^\k\ + ¥) + (7 

(64) 

Here A; is a quantum number of relativistic angular momentum: k = —I, l + l, I is a quantum 
number of orbital moment. Comparison of right-hand sides of formulas (163]) and fl64l) yields 
that AMqo = lim^oo AM n (0) = m is a rest mass m of an electron. Hence, Eq. fl6"31 is an 
equation of motion for a relativistic electron - hole pair with a reduced mass = m/2 

which, apparently, is a relativistic charged vector-boson considered in quasirelativistic limit 
n — > oo. 



6 Conclusion 

The method of Hartree - Fock self-consistent field in non-relativistic case or of Dirac - Fock 
self-consistent field in relativistic case is applicable for description of valent electron and 
polarization of atomic core only within the framework of the perturbation theory with the 
pseudopotential method as a zero-order approximation. To describe effects of polarization 
of electron density distribution in the case of intensive electromagnetic interactions with a 
atomic core it is necessary to take into account correlation interaction which is neglected in 
the above mentioned method. The method offered in this paper allows to describe effects of 
polarization not only for a atomic core in a field of valent electron, but also polarization of 
atom as the whole in the electron-hole formalism. 
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